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Abstract 

We investigated ram-pressure effects by an 
intracluster wind on an inner disk of spi- 
ral galaxies by a hydrodynamical simulation. 
Even if the wind is mild and not strong enough 
to strip the gas disk, the ram pressure disturbs 
orbits of the inter-arm gas significantly. This 
results in asymmetric dense molecular arms in 
the inner few kpc region of a galaxy. This 
mechanism would explain the asymmetric CO 
gas distributions in the central regions often 
observed in Virgo spirals. 

Key words: galaxies: spiral — galaxies: 
kinematics and dynamics — galaxies: ISM — 
galaxies: cluster of — intergalactic matter 

1 Introduction 

Ram-pressure by intracluster medium (ICM) 
causes the stripping of interstellar matter 
(ISM) from galaxies (Farouki, Shapiro 1980; 
Kritsuk 1984; Gaetz et al. 1987; Balsara et 
al. 1994; Sofue 1994; Quilis et al. 2000). It 
also produces a disturbed distribution of ISM 
in galaxies, such as head-tail H i structures 
as observed in Virgo galaxies (Cayattc et al. 
1990; VoUmer et al. 2000, 2001; Phookun et 
al. 1993, 1995). The ram-pressure effect has, 



thus, been discussed mainly in relation to H i 
gas stripping and outer disk structures. How- 
ever, little attention has been paid to its effect 
on the inner molecular disk and arms. Only a 
few authors have discussed the ram effect on 
the spiral structure (Tosa 1994) and molecular 
clouds (Kenney et al. 1990; Sofue 1994). 

Kenney et al. (1990) have shown that Virgo 
galaxies often exhibit asymmetric inner molec- 
ular disks. A recent high- resolution CO-line 
survey of Virgo cluster galaxies such as NGC 
4254 and NGC 4654 (Sofue et al. in prepara- 
tion) has revealed that some of them show a 
significant asymmetry of the inner molecular 
disks and arms. It is not clear if such an in- 
ner deformation of dense gas disks can be pro- 
duced by ram pressure effect, which is thought 
to be the cause of their deformed H i envelops. 
Since these two galaxies have no massive com- 
panion that can disturb such inner disks, and 
they both show a head-tail H i outer structure 
(Phookun et al. 1993, 1995), the inner defor- 
mation of molecular disks could also be due 
to the ram-pressure effect, while such an inner 
ram effect has not yetbeen investigated. 

If the wind is very strong, such as assumed 
by Quilis et al. (2000) for the core of a rich 
cluster with an ICM density of ~ 3 x 10~^ 
atoms cm^'^ and a wind velocity higher than 
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~ 2000 km , the ISM of any galaxies would 
be completely stripped. On the other hand, if 
the wind is mild, in such a case as for the Virgo 
cluster, where the ICM density is of the order 
of IQ-^ - IQ-^ and the velocity is ~ 1000 km 
s"-'^ , outer H i envelopes are deformed to pro- 
duce head-tail structures (VoUmer et al. 2000, 
2001). In current simulations, such as those by 
Vollmcr et al., which were aimed at gas strip- 
ping and tailing of the outer H i disks, the 
detailed behavior of the inner disk gas inside 

10 kpc was not well understood because of 
the resolution. 

In the present paper, we consider a mild 
ICM wind, and discuss its effect on the inner 
disk gas based on 2D hydrodynamical simula- 
tions with higher resolution than those aimed 
at outer H i stripping, as above. If we sim- 
ply apply the ram-stripping condition to an 
azimuthally structure-less gas disk, the ram 
pressure would hardly affect the inner disk. 
However, if we consider a spiral structure with 
an arm-to-interarm density contrast, it may 
happen that the ram-prcssurc can affect the 
low-density interarm gas. We consider here 
the possibility that the ram-pressure can af- 
fect the dense molecular gas within the central 
few kpc region of spiral galaxies through dis- 
turbances of the orbits of inter-arm gas, even 
if the ICM wind is not strong enough to strip 
the disk gas. 



2 Ram-Pressure Force on 
the Arm and Interarm 
Gases 

The component of the ram force parallel to the 
galactic plane exerted by an intergalactic wind 
on a gas element is given by 

/iCM ^ PiCM s Sv^ cosa sina, (1) 



where picM and pisM are the gas densities 
of ICM and ISM, respectively, s is the sur- 
face area of the element, and a is the angle 
between the wind direction and the galactic 
plane (Farouki, Shapiro 1980; Kritsuk 1984; 
Sofue 1994; Tosa 1994). The motion of the 
undisturbed ISM element is governed by the 
gravitational force, which is approximately 
equal to the centrifugal force, 

/ism ~ PiSM d s v'^^^/R, (2) 

where d is the thickness of the gas disk, and 
R is the galactocentric radius of the clement. 
Now, the ratio of /icm to /ism is given by 

nicM R f dv .2 . 
r] 1 1 cosa sma, (3) 

niSM d ^Vrot 

where nicM and tiism are the number densi- 
ties of hydrogen atoms of the ICM and ISM, 
respectively. If rj exceeds unity, the ram force 
can disturb the ISM motion, while if it is 
smaller than unity, the ISM motion is little 
affected. 

Let us consider the inner part of a galaxy 

at i? ~ 5 kpc with d ^ 100 pc, rotating at 
Viot ~ 200 km s~^ . For an ICM wind with 
?^iCM ^ 10~^ cm~^, VicM ~ 10^ km , and 
a 45°, we obtain the ratio to be 

V ~ 0.6 nfs^M [H cm-3]. (4) 

This relation implies that the ISM is stripped 
if nisM 1 H cm^"*, since the force perpendic- 
ular to the galactic plane is of the same order. 
We stress, however, that the relation indicates 
that the orbits of the gas within the disk plane 
is significantly disturbed if njsM ~ 1 H cm""^, 
even if the wind is not strong enough to strip 
the gas. This may indeed apply to the inter- 
arm ISM in the inner disk within a few kpc 
radius. On the other hand, high-density galac- 
tic shocked arms, where nisM ^ 1 H cm~^, 
would be hardly disturbed. We, thus, antici- 
pate that, even if the ICM wind is mild, not 
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strong enough to strip the disk, the interarm 
gas in the inner disk would be significantly dis- 
turbed, which may result in deformed galactic 
shock waves. Since 6v is greater on the head- 
wind side of the rotation axis compared to the 
following-wind side, the deformed shock waves 
could be asymmetric with respect to the rota- 
tion axis. 

Figure 1 illustrates the ram-deformation 
mechanism of the dense molecular arms, and 
how a lopsided spiral pattern is created. Since 
the density of the inter-arm gas (A and A' in 
Fig. 2) is much lower than the average den- 
sity of the disk, the ram force by the ICM wind 
(thin lines) easily disturbs the orbits of inter- 
arm gas (thick curved arrows). The gas on 
the distorted orbits encounters density waves 
(dashed spirals) at different places (B and B') 
from those expected for undisturbed orbits 
(dashed lines), and produces deformed dense 
molecular arms (thick spirals). In the next 
sections we discuss a numerical simulation of 
the ram deformation of spiral arms in order 
to understand whether this mechanism can in- 
deed create deformed shocked arms, and how 
the deformed arms look like in realistic model 
disks. 

— Fig. 1 — 

3 Numerical Simulation 

3.1 Bcisic assumptions 

The ram-pressure acceleration per unit mass 
is given by 



Cn 



ICM 



ViCM — Vrot 



Viot 



)■ 

(5) 



where nicM is the ICM density, Vicm is the 
ICM velocity with respect to the galaxy, Vyot 
is the circular velocity of the clouds, and C is 
evaluated to be on the order of (e.g., Tosa 
1994). 



We assume that the galactic disk is thin 
and faces the ICM wind everywhere, being 
not shielded by neighboring clouds. Because 
this assumption would not apply to an edge- 
on wind, we consider here a wind with a mild 
inclination. We consider the inner disk of a 
galaxy, where the gravitational potential is 
deep and the ISM is dense enough so that 
stripping does not occur, as discussed in the 
previous section, and we treat a 2D disk in a 
fixed rotating potential. 

For the velocity of a galaxy in the ICM Vicm 
and the ICM density nicM, we adopt three val- 
ues for each parameter. The galaxy's veloc- 
ity Vicm, is taken to be 530, 1000, and 1500 
km s^^ , where the first value is suggested by 
Phookun, Mundy (1995) for NGC 4654 in the 
Virgo cluster. The ICM density, Uicm is taken 
to be 1 X 10-* cm-3, 5 x lO"* cm'^, and 1 x IQ-^ 
cm~^, where the first value is typical for the 
intergalactic density. 

3.2 Numerical Method 

For simplicity, we assumed that the interstel- 
lar gas is ideal, inviscid, and compressible. 
We used a freely downloadable and usable 
hydrodynamical code, VH-1 (Blondin, Lufkin 
1993). This is a multidimensional hydrody- 
namics code for an ideal compressible fluid 
written in FORTRAN, developed by the nu- 
merical astrophysics group at the University 
of Virginia based on the Picccwise Parabolic 
Method Lagrangian Remap (PPMLR) scheme 
of Colella, Woodward (1984). The PPMLR 
has the advantage of maintaining contact dis- 
continuities without the aid of a contact steep- 
ener, and is sufficiently good to be applied 
to a galactic-scale hydrodynamical simulation. 
The code does not take into account the 
gas's self-gravity, artificial viscosity, variable 
gamma equation of state, and radiative heat- 
ing and/or cooling. The self-gravity of the gas 
is not taken into account, because we consider 
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a case where the gas mass is not so much as to 
contribute to the density wave potential by the 
stellar disk, and also because we do not intend 
to discuss such processes as clumping of gas 
and cloud formation in the arm and inter- arm 
regions. The interstellar gas can be assumed 
to be isothermal, as is assumed here, while if 
cooling is taken into account, the shocked gas 
arms would become much denser than those 
calculated below. However, all such neglec- 
tion would not affect the physical essence of 
the present study, which was aimed at simu- 
lating ram-deformation of the orbits of inter- 
arm gas and resulting disturbed shocked arms, 
as illustrated in figure 1. 

3.3 Gravitational Potential 

We implicitly give a gravitational potential, 
which comprises the following two terms: (i) a 
static axisymmetric potential, and (ii) a non- 
axisymmetric, rotating bar potential. The po- 
tential is expressed by 

^R,,j>) = MR) + ^i{RA)- (6) 

We adopt a "Toomre disk" (Toomre 1981) 
potential for the axisymmetric component as 
given by 



MR) = - 



1 



a {R? + o2)V2 ' 



(7) 



where a is the core radius and c = 
(27/4)i/4a. Through our numerical sim- 
ulation, we fixd the core radius and maxi- 
mum circular velocity to be a = ^/2 kpc and 
t'max = 200 km s~^ . 

The nonaxisymmetric potential was taken 
from Sanders (1977), assuming rigid rotation 
at a pattern speed, rip, which has the form 



(i?2 + a2)3/2 



$o(-R)cos2(0-Oj,i) 
(8) 



where e is the strength of the bar of the order 
of £ = 0.15. Spiral shocked arms of gas are 
produced by this potential. 

3.4 Initial Conditions 

Initially, we set 256 x 256 two-dimensional cells 
corresponding to 12.8 kpc xl2.8 kpc field, 
while setting the field center at the coordinates 
origin. The initial number density was taken 
to be 5 cm~^ the inner disk at i? > 8 kpc disk, 
and 1 cm"-'^ at i? > 8 kpc. The initial rota- 
tion velocity of each gas cell was set so that 
the centrifugal force would balance the gravi- 
tation. The bar pattern speed, fip, was taken 
to be 23 km s~^ , and the strength of the bar 
e was taken to be 0.10. 



4 Ram-Pressure Deforma- 
tion of Dense Molecular 
Arms 

4.1 Deformation of Inner Spiral 
Structure 

Figure 1 shows the result of a simulation in- 
cluding ram-pressure effects for the various pa- 
rameter combinations, as described in the pre- 
vious section. Both the ISM and spiral pattern 
rotate counterclockwise, and the ICM wind 
blows from left to right. The simulation shows 
that the orbits of the diffuse inter-arm gas are 
easily disturbed by the ram force, which re- 
sults in a significant displacement of galactic 
shock waves from their undisturbed symmetric 
positions, as illustrated in figure 1. 

Highly asymmetric dense spiral arms in the 
central region are produced by this mechanism 
if the wind speed is higher than ~ 1000 km 
s~^ and the ICM density is greater than sev- 
'eral 10~^ H cm~^. A head-tail structure of 
dense gases slanted to the ICM wind, like NGC 
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4654 nucleus, can be produced by this mecha- 
nism, if nicM X V^cM is greater than ~ 3 x 10^^ 
cm~^ s^^. One spiral arm on the downstream 
side is prominent, reproducing the lopsided 
arms, as observed in NGC 4254 and NGC 
4654. 

— Fig. 2 — 

4.2 Deformed Molecular Arms 

We then calculated the distribution of molec- 
ular fraction (Elmegreen 1993) corresponding 
to figure 1. The molecular fraction is defined 
by 

PH2 2nH2 /ri\ 

/mol = ; = ; , (9j 

Phi + Ph2 nm + 

where pu2, Pm, and tihi are the mass 
and number densities of molecular and H i 
gases, respectively. We used a method de- 
scribed by Sofue et al. (1995) and Honma et al. 
(1995) to calculate the galaxy-scale molecular 
fraction; they investigated molecular fronts in 
spiral galaxies using a phase-transition model 
proposed by Elmegreen (1993). In this model, 
the molecular fraction is determined by three 
parameters; the interstellar pressure, P, the 
UV radiation field, U, and the metallicity, Z. 
For U and Z, we adopted an exponential func- 
tion of galactocentric radius, and calculated 
/mol for corresponding gas pressure, P, which 
is determined by the gas density in each cell. 

Figure 3 shows the result of a numerical 
simulation of the molecular fraction. The in- 
ner few kpc region is dominated by molecular 
gas, where the molecular fraction is as large 
as ~ 70 — 90%. Also, the molecular frac- 
tion increases suddenly at the galactic shocks, 
which are already deformed from symmetric 
arms. Thus, the simulation has revealed that 
highly deformed inner molecular arms can be 
produced by ram-pressure disturbances on the 
inter-arm low-density regions. 

— Fig. 3 — 



4.3 Comparison with Observa- 
tions and Wind Velocities 

We now compare the results with H I and CO- 
line observations of the Virgo galaxies, NGC 
4254 and NGC 4654. The heliocentric radial 
velocity of NGC 4254, 2407 km s'^ , is about 
1100 kms-i different from that of NGC 4486 
(M 87), the center of the Virgo cluster, 1282 
km s~^ (de Vaucouleurs et al. 1991). Also, by 
comparing with the lopsided H i distribution 
of NGC 4254 (Phookun et al. 1993) and the 
ram pressure simulation on galaxies (Abadi et 
al. 1999), we can exclude the possibility of any 
face-on motion of NGC 4254. Assuming that 
NGC 4254's orbit is inclined by 45° from the 
line of sight, we may estimate the velocity of 
motion of NGC 4254 in the Virgo cluster as 
being ~ 1500 km s~^ . Then, the rotation 
velocity of NGC 4254 (~ 150 km s'^ ) is neg- 
ligible compared to the wind velocity, so that 
the H I tail grows toward downstream. The 
location of a prominent spiral arm (lye et al. 
1982) and the direction of rotation are consis- 
tent with our numerical simulation. However, 
if we assume a slower wind velocity, e.g., on 
the order of, or smaller than, ~ 750 km s~^ , 
the simulation cannot reproduce the observed 
features. 

Our simulation for moderate ICM veloc- 
ity, which predicts an off-center bar of dense 
gas tilted toward the ICM wind, is consis- 
tent with the observations of NGC 4654 in 
H I (Phookun, Mundy 1995) and CO (Sofue 
et al. in preparation). Our simulations show 
that spiral arm in the upstream side becomes 
stronger than that in the downstream side. 
Considering the location of the prominent op- 
tical arm (Frei et al. 1996) and an elongation 
of the observed molecular bar and H i tail, we 
find that the direction of motion of NGC 4654 
is toward the northwest, with a significantly 
high velocity compared to the velocity disper- 
sion of the Virgo cluster. Taking it into ac- 



5 



count that the rotation velocity of NGC 4654 
is relatively slow, the velocity of NGC 4654 in 
the Virgo cluster would be greater than 1000 
km s~^ , although the heliocentric radial ve- 
locity, 1054 km s~^ , is close to Virgo's central 
velocity. 

5 Discussion 

We have considered the ram-pressure effects of 
a mild ICM wind on gaseous disks in cluster 
galaxies. Galaxies in the central region of the 
Virgo galaxies show H i deficiency, where the 
molecular gas fraction is higher than that of 
the galaxies in the outer region of the cluster 
(Kenncy ct al. 1990). This fact has been nat- 
urally explained by the ram-prcssTirc stripping 
of the H I gas; the ram effects arc negligible 
on the inner molecular disks, while it is cru- 
cial for the H i outer disk and envelops. On the 
other hand, it is also known that ram-affected 
Virgo galaxies show asymmetric CO gas dis- 
tributions (Kenney et al. 1990; Sofue et al. in 
preparation). 

In the present paper, we have shown that 
the orbits of the inter-arm ISM are disturbed 
by the ram pressure of the ICM wind, even 
if the wind is mild and not strong enough to 
strip the gas disk. The disturbed inter-arm 
ISM causes highly asymmetric moleciilar arms 
in the inner few kpc of the disk. The 3D sim- 
ulation by VoUmer et al. (2000, 2001) for 
a similar wind condition shows that the in- 
ner disk within 10 kpc radius is not stripped, 
but suffers from perpendicular disturbances to 
the disk plane. Vertical displacements may 
result in off-plane molecular structure in the 
inner few kpc disk, in addition to the asym- 
metric arms as simulated here using the 2D 
scheme. Also, such 3D effects as the Kelvin- 
Helmholtz instability by the shearing motion 
between the disk and ICM (e.g. Mori and 
Burkert 2000) cannot be touched upon by the 



present 2D simulation. Therefore, detailed 3D 
simulations will be crucial to thoroughly un- 
derstand the inner ram effect in more detail, 
while the present 2D results can tell us about 
some essential mechanism to cause the non- 
axisymmetric molecular structures in the in- 
ner disk. 

However, in the cases of much stronger 

winds, the 2D assumption cannot be applied 
in any way, and 3D treatment of stripping is 
crucial. In fact, ram effect by a wind with a 
pressure much higher than that considered in 
this paper has been simulated by Quilis et al. 
(2000) with a 3D hydro dynamical simulation; 
they have shown that the ISM in a galaxy is 
completely stripped within one galactic rota- 
tion. 
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orbits (thick arrows) encounters density waves 
(dashed spirals) at different places (B and B') 
from those expected for undisturbed orbits 
(dashed spirals), and produces deformed dense 
gaseous arms (thick spirals). 

Fig. 2. Snapshots of the density distri- 
bution in the ram-pressure models after 1.3- 
1.7 X 10* yr. The labels indicated the adopted 
parameters. Alphabets: ICM velocity, Vicm = 
(a) 530, (b) 1000, and (c) 1500 km s'^ . 
Roman numbers: ICM density njcM = (i) 
1 X 10"'', (ii) 5 X 10"'', and (iii) 1 x lO'^ cm'^. 
The color-density key is shown at the bottom, 
which is common to all snapshots. The rota- 
tion direction of gases and a bar potential is 
counterclockwise, and ICM wind brows from 
left to right. The strength of the bar e is 0.10 
in all models. Weak straight waves toward up- 
per and right boundaries in b-(ii, iii) and c-(ii, 
iii) arc artifact due to numerical reflection at 
the boundaries 

Fig. 3. Same as figure 2, but showing the 
distribution of molecular fraction, /moi- 



Figure Captions 

Fig. 1. Schematic illustration of the ram- 
deformation mechanism to cause asymmetric 
inner molecular arms. Orbits (thin dashed ar- 
rows) of low-density inter-arm gas (A and A') 
are disturbed by the ram pressure of the ICM 
wind (thin arrows). The gas on the distorted 
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